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<N 

O ' We consider the time evolution of a system of N identical bosons whose interaction 

potential is rescaled by TV -1 . We choose the initial wave function to describe a condensate 
(-^ ' in which all particles arc in the same one-particle state. It is well known that in the mean- 

field limit N — * oo the quantum iV-body dynamics is governed by the nonlinear Hartree 
equation. Using a nonperturbative method, we extend previous results on the mean-field 
limit in two directions. First, we allow a large class of singular interaction potentials as well 
' as strong, possibly time-dependent external potentials. Second, we derive bounds on the 

rate of convergence of the quantum iV-body dynamics to the Hartree dynamics. 

J: 
+-> 

1. Introduction 

. We consider a system of iV identical bosons in d dimensions, described by a wave function 

m ! * N € n {N) . Here 

# := Ll(R Nd ,d Xl ...dx N ) 

is the subspace of L 2 (M Nd , dx\ ■ ■ • dx/v) consisting of wave functions ^>n(xi, . . . ,xn) that are 
symmetric under permutation of their arguments x±, . . . , xn G R rf - The Hamiltonian is given by 

ON ! N 



O 



Hn = J2 h i + Jj ^2 w ( x i~ x j), ( L1 ) 

where hi denotes a one-particle Hamiltonian h (to be specified later) acting on the coordinate 
Xi, and w is an interaction potential. Note the mean- field scaling 1/N in front of the interaction 
potential, which ensures that the free and interacting parts of Hn are of the same order. 
The time evolution of is governed by the iV-body Schrodinger equation 

i9ttfjv(t) = H N * N (t), *at(0) = *tv,o- (1.2) 



For definiteness, let us consider factorized initial data ^n,o = f° r some G L (R ) satisfy- 
ing the normalization condition ||¥>o||z, 2 (R d ) = 1- Clearly, because of the interaction between the 
particles, the factorization of the wave function is not preserved by the time evolution. However, 
it turns out that for large N the interaction potential experienced by any single particle may be 
approximated by an effective mean- field potential, so that the wave function *$>N(t) remains ap- 
proximately factorized for all times. In other words we have that, in a sense to be made precise, 
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vPjv(i) ~ for some appropriate (p(t). A simple argument shows that in a product state 

tp(t)® N the interaction potential experienced by a particle is approximately w * \ (p(t)\ 2 , where * 
denotes convolution. This implies that <p{t) is a solution of the nonlinear Hartree equation 

idMt) = fup(t) + (w*\<p(t)\ 2 )<p(t), <p(0) = <p . (1.3) 

Let us be a little more precise about what one means with ~ cp® N (we omit the irrelevant 
time argument). One does not expect the L 2 -distance \\*&n — ( P^ N \\ L 2^Nd^ to become small as 
N — > oo. A more useful, weaker, indicator of convergence should depend only on a finite, fixecfl 
number, k, of particles. To this end we define the reduced fc-particle density matrix 

7 j? := Ttfc+i jv l^/vX^/vl , 

where Trfc+i,,,,^ denotes the partial trace over the coordinates x k+ i, . . . ,xjy, and l^jv^Vl/jvl 
denotes (in accordance with the usual Dirac notation) the orthogonal projector onto \&at. In 
other words, 7 jy is the positive trace class operator on L^(M kd , dxi ■ ■ ■ dx k ) with operator kernel 

jP(x 1 ,...,x k ;y 1 ,...,y k ) = / dx k+1 ■ ■ ■ dx N *jv(xi, • • • , x N )^ N (yi, . . . , y k , x k+1 , . . . , x N ) . 



(k) 

The reduced £>particle density matrix 7^ embodies all the information contained in the full 
iV-particle wave function that pertains to at most k particles. There are two commonly used 
indicators of the closeness yffi « (\ip) ((p\)® k : the projection 

(fc) 1 /,„®fe A k ),„®k\ 



and the trace norm distance 



:= Tr 



7^-(M(y|) M • (i-4) 



It is well known (see e.g. [9]) that all of these indicators are equivalent in the sense that the 
vanishing of either jR$ or for some k in the limit N — > 00 implies that lim^v ^ = 
lim^v E N = for all k . However, the rate of convergence may differ from one indicator to 
another. Thus, when studying rates of convergence, they are not equivalent (see Section [2] below 
for a full discussion). 

The study of the convergence of Jpj\t) in the mean-field limit towards (\tp(t)} ((p(t)\)® k for 
all t has a history going back almost thirty years. The first result is due to Spohn [13], who 



showed that limiv R^\t) = for all t provided that w is bounded. His method is based on the 
BBGKY hierarchy, 



k 



N 

i=l l^i<j^fc 

+ N N k E Tffc+i [w(xi - x k+1 ),^ +1 \t)] , (1.5) 

i=l 



1 In fact, as shown in Corollary 13.21 k may be taken to grow like o(N). 
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an equation of motion for the family (7L (t))fceN of reduced density matrices. It is a simple 
computation to check that the BBGKY hierarchy is equivalent to the Schrodinger equation (II. 2p 
for \I/jv(t). Using a perturbative expansion of the BBGKY hierarchy, Spohn showed that in the 
limit N — > 00 the family (7^^(i))fceN converges to a family (7S (i))feeN that satisfies the limiting 
BBGKY obtained by formally setting N = 00 in (jl.5p . This limiting hierarchy is easily seen to 
be equivalent to the Hartree equation (j 1 .3f) via the identification 7^»' ) (i) = (\ip(t))(ip(t)\)^ k . We 
refer to [3] for a short discussion of some subsequent developments. 

In the past few years considerable progress has been made in strengthening such results in 
mainly two directions. First, the convergence lirn/v-R^r (i) = for all t has been proven for 
singular interaction potentials w. It is for instance of special physical interest to understand the 
case of a Coulomb potential, w(x) = A|x| _1 where A £ R. The proofs for singular interaction 
potentials are considerably more involved than for bounded interaction potentials. The first 
result for the case h = —A and w(x) = A|x| _1 is due to Erdos and Yau [3]. Their proof uses the 
BBGKY hierarchy and a weak compactness argument. In pQ, Schlein and Elgart extended this 
result to the technically more demanding case of a semirelativistic kinetic energy, h = VI - A 
and w{x) = A|x| _1 . This is a critical case in the sense that the kinetic energy has the same 
scaling behaviour as the Coulomb potential energy, thus requiring quite refined estimates. A 
different approach, based on operator methods, was developed by Frohlich et al. in [3], where 
the authors treat the case h = —A and w(x) = A|x| _1 . Their proof relies on dispersive estimates 
and counting of Feynman graphs. Yet another approach was adopted by Rodnianski and Schlein 
in p2]- Using methods inspired by a semiclassical argument of Hepp [6] focusing on the dynamics 
of coherent states in Fock space, they show convergence to the mean-field limit in the case 
h = —A and w(x) = A|x| _1 . 

The second area of recent progress in understanding the mean-field limit is deriving estimates 
on the rate of convergence to the mean-field limit. Methods based on expansions, as used in |13] 
and [3J, give very weak bounds on the error (t), while weak compactness arguments, as used 
in [3] and [JJ, yield no information on the rate of convergence. From a physical point of view, 
where N is large but finite, it is of some interest to have tight error bounds in order to be able 
to address the question whether the mean-field approximation may be regarded as valid. The 
first reasonable estimates on the error were derived for the case h = —A and w(x) = \\x\~ 1 
by Rodnianski and Schlein in their work [12] mentioned above. In fact they derive an explicit 
estimate on the error of the form 



for some constants C\(k), C2(k) > 0. Using a novel approach inspired by Lieb- Robinson bounds, 
Erdos and Schlein [2] further improved this estimate under the more restrictive assumption that 
w is bounded and its Fourier transform integrable. Their result is 



for some constants C\, C2, C3 > 0. 

In the present article we adopt yet another approach based on a method of Pickl [10j . We 
strengthen and generalize many of the results listed above, by treating more singular interaction 
potentials as well as deriving estimates on the rate of convergence. Moreover, our approach 




Cl ( k ) c C 2 (k)t 



3 



allows for a large class of (possibly time-dependent) external potentials, which might for instance 
describe a trap confining the particles to a small volume. We also show that if the solution ip{-) 
of the Hartree equation satisfies a scattering condition, all of the error estimates are uniform in 
time. 

The outline of the article is as follows. Section [2] is devoted to a short discussion of the 
indicators of convergence E$ and , in which we derive estimates relating them to each 
other. In Section [3] we state and prove our first main result, which concerns the mean-field limit 
in the case of L 2 -type singularities in w; see Theorem 13.11 and Corollary 13.21 In Section U] we 
state and prove our second main result, which allows for a larger class of singularities such as 
the nonrelativistic critical case h = —A and w(x) = \\x\~ 2 ; see Theorem 14.11 For an outline of 
the methods underlying our proofs, see the beginnings of Sections [3] and HI 

Acknowledgements. We would like to thank J. Frohlich and E. Lenzmann for helpful and stim- 
ulating discussions. We also gratefully acknowledge discussions with A. Michelangeli which led 
to Lemma 12. 11 

Notations. Except in definitions, in statements of results and where confusion is possible, we 
refrain from indicating the explicit dependence of a quantity a^(t) on the time t and the particle 
number N. When needed, we use the notations a(t) and a\t interchangeably to denote the value 
of the quantity a at time t. The symbol C is reserved for a generic positive constant that may 
depend on some fixed parameters. We abbreviate a ^ Cb with a < b. To simplify notation, we 
assume that t ^ 0. 

We abbreviate L p (M. d ,dx) = LP and ||-||lj> = ||-|| p . We also set || • || J L 2 (R JVd ) = II "II- For set 
we use H s = H s (R d ) to denote the Sobolev space with norm ||/||_fp = ||(1 + |^| 2 ) s ^ 2 /|| 2 , where 
/ is the Fourier transform of /. 

Integer indices on operators denote particle number: A A;-particle operator A (i.e. an operator 
on H^) acting on the coordinates Xi l , . . . ,Xi k , where i\ < ■ ■ ■ < is denoted by A^.,,^. Also, 
by a slight abuse of notation, we identify /c-particle functions f(x\, . . . , Xk) with their associated 
multiplication operators on Ti^. The operator norm of the multiplication operator / is equal 
to, and will always be denoted by, ||/||oo- 

We use the symbol Q(-) to denote the form domain of a semibounded operator. We denote 
the space of bounded linear maps from X\ to X2 by C{X\\ X2), and abbreviate C(X) = C(X; X). 
We abbreviate the operator norm of C(L 2 (R Nd )) by 1 1 - 1 1 . For two Banach spaces, X\ and X2, 
contained in some larger space, we set 



This section is devoted to a discussion, which might also be of independent interest, of quanti- 
fy) (k) 

tative relationships between the indicators E N and R N . Throughout this section we suppress 




H/ll-XinXa 



inf (||A|| Xl + ||/ 2 || X2 ) 
t=h+h 

ll/lk + ||/||x 2) 



and denote by X\ + X2 and X\ n X2 the corresponding Banach spaces. 



2. Indicators of convergence 
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the irrelevant index N. 

Take a A;-particle density matrix 7W G C(7i^) and a one-particle condensate wave function 
</? G L 2 . The following lemma gives the relationship between different elements of the sequence 
EW,E( 2 \..., where, we recall, 

£?(*) = 1-^8*,/)^). (2.1) 

Lemma 2.1. lei 7^) G safe/y 

7 {fc) > 0, Tr 7 (fe) = 1. 
Let (p G L 2 satisfy ||y>|| = 1. TTien 

< kE^K (2.2) 
Proof. Let (^ , f ) ) i>1 be an orthonormal basis of H^ k) with <f>[ k] = Then 

i^l i>2 

= , 7 (1) ^>-£<^ _1) , 7 (fe) ¥> ® ~ 1} > • 

i>2 

Therefore, 

(^,7 (1) ^)-<^ fc ,7 (fc) ^ fc > 
= ^< V ®<6i*- 1) ,7 ()fc) V®* i ( *~ 1) ) 

i>2 j>l 

= E E ® *< (fr_1) > *f ® *? -1) > - E (*f ® ^ (fc_1) > *? } ® 

= l_^®(fc-l) )7 (fc-l)^®(fc-l)). 

This yields 

and the claim follows. □ 

Remark 2.2. The bound in (|2.2p is sharp. Indeed, let us suppose that E^ <c k f(k) E^ for 
some function /. Then 

E^ l-(l-a) k l-fl-al* 
f(k) 2 sup > sup ^ ^ > lim \ >- = 1, 

7 ( fe ) kEW o<a<i ka a-*o ka 

where the second inequality follows by restricting the supremum to product states 7W = 
and writing a = E^ l \ 
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The next lemma describes the relationship between E^ and R^ k \ where, we recall, 

R (k) = it _ (\<p)(<p\ f k \ . 

Lemma 2.3. Let j( k > £ C(H^) be a density matrix and (p £ L 2 satisfy \\ip\\ = 1. Then 

E^ < (2.3a) 
i?( fe ) ^ y/8EW . (2.3b) 

Proof. It is convenient to introduce the shorthand 

p(*> := (\<p){<p\)° k . 

Thus, 

# = l-(/ )7 W^) = Tr(p( fc )-pW 7 (fc) ) < b (fe) ||TV|pW- 7 W| = 

which is (|2.3aj) . In order to prove (|2.3b|) it is easiest to use the identity 

Tr |p( fe ) -7 (fc) | = 2 — , (2-4) 

valid for any one-dimensional projector p^ and nonnegative density matrix 7W. This was first 
observed by Seiringer; see [12]. For the convenience of the reader we recall the proof of (|2.4p . 
Let (A n ) ne pj be the sequence of eigenvalues of the trace class operator A := 7W — p( k \ Since 
is a rank one projection, ^4 has at most one negative eigenvalue, say Ao- Also, Tr A = 
implies that J2 n ^« = ^- Thus, X^nl-^l = 2 1 Ao | , which is (|2.4j) . 
Now ([Z3D yields 

flW = Tir|p( fc )- 7 (*)| = 2||p( fc )- 7 ( fc )|| < 2^Tr(pW - 7 W) 2 . 

Then (I2.3b[) follows from 

Tr (p( fc )- 7 ( fc )) 2 = l-2Tr(p( fc ) 7 W) +Tr( 7 «) 2 < - Tr(p^^) + 1 = 2E™ . 

Alternatively, one may prove (|2~3b|) without (p^|) by using the polar decomposition and the 
Cauchy-Schwarz inequality for Hilbert-Schmidt operators. □ 



Remark 2.4. Up to constant factors the bounds (12. 3p are sharp, as the following examples show. 
Here we drop the irrelevant index k. Consider first 

1\ fl-a 



* = ' 7 V a 



where ^ a ^ 1. As above we set p := \(p){<p\. One finds 

E = 1 — (tp^tp) = a, R = Tr\p — 7 | = 2a. 
so that (I2.3ah is sharp up to a constant factor. 
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It is not hard to see that if 7 and p commute then (|2.3bj) can be replaced with the stronger 
bound R < E. In order to show that in general (I2.3b[) is sharp up to a constant factor, consider 

l\ ( 1 — a V« — a? 

2 



V' 7 - \Va- 

where ^ a ^ L One readily sees that 7 is a density matrix (in fact, a one-dimensional 
projector). A short calculation yields 

E = 1 - (<p,7<p) = a 

as well as 

Tr (7(1 -p)\ = y/a. 

Using 

TV 1 7(1— p) I = Tr|7 — p + p — 7p| ^ 2 Tr|p — 7I 

we therefore find 



as desired. 



m 1 \/i VE 

R = Trp-7 > — = — 
2 2 



3. Convergence for L 2 -type singularities 

This section is devoted to the case w S L 2 + L°° . 

3.1. Outline and main result. Our method relies on controlling the quantity 

a N (t) := E$(t). (3.1) 
To this end, we derive an estimate of the form 

d N (t) ^ A N (t) + B N (t)a N (t), (3.2) 
which, by Gronwall's lemma, implies 

t ft t 

a N (t) ^ a N (0)efo B N+ / A N (s) Bn ds . (3.3) 

^0 

In order to show (|3,2p . we differentiate ajv(i) and note that all terms arising from the one- 
particle Hamiltonian vanish. We control the remaining terms by introducing the time-dependent 
orthogonal projections 

pit) := \<p(t))(<p(t)\, q(t) := 1 - p(t) . 

We then partition 1 = p(t)+q(t) appropriately and use the following heuristics for controlling the 
terms that arise in this manner. Factors p(t) are used to control singularities of w by exploiting 
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the smoothness of the Hartree wave function <p(t). Factors q(t) are expected to yield something 
small, i.e. proportional to ajv(t), in accordance with the identity ajv(i) = (\l/jv(£) , (/i(£)\Pjv(£))- 
For the following it is convenient to rewrite the Hamiltonian as 



N 



H N = J2h i + ± W n =■ H°N + H% , (3.4) 



N 



where Wij := w(xi — xj). We may now list our assumptions. 

(Al) The one-particle Hamiltonian h is self-adjoint and bounded from below. Without loss of 
generality we assume that h ^ 0. We define the Hilbert space X N = Q(H%) as the form 



domain of with norm 



1*11^ == ||(l + #ft) 1/2 *||. 



(A2) The Hamiltonian (|3.4f) is self-adjoint and bounded from below. We also assume that 
Q(H N ) c X N . 

(A3) The interaction potential w is a real and even function satisfying w £ L Pl + L P2 , where 
2 ^ Pi ^ P2 ^ °°- 

(A4) The solution (p(-) of (|1.3p satisfies 

where 2 ^ (72 ^ <Zi ^ oo are defined through 

111 , . 

- = - + -, i = 1,2. 3.5 

2 ^ ft 

Here X]" denotes the dual space of X\, i.e. the closure of L 2 under the norm := 

We now state our main result. 

Theorem 3.1. Let ^/ N:0 e Q(H N ) satisfy \\^n,o\\ = 1? and tp G X 1 n L qi satisfy \\ip \\ = 1. 
Assume that Assumptions (Al) - (A4) /ioZ<i. TTien 

4>{t) := 32\\w\\ LPl+LP2 [ ds(\\cp(s)\\ gi + \\tp(s)\\ q2 ). 

J 



where 



'0 

We may combine this result with the observations of Section [2j 
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Corollary 3.2. Let the sequence ^?n,o G Q{Hn), N G N, satisfy the assumptions of Theorem 
I3.il as uieZ/ as 

Taen we have 

Ef{t) < A e «*>, < y^e^)/ 2 . 

Remark 3.3. Corollary 13.21 implies that we can control the condensation of k = o(N) particles. 

Remark 3.4. Assumption (A3) allows for singularities in w up to, but not including, the type 
| ^ | —3/2 - n ^jy-gg dimensions. In the next section we treat a larger class of interaction potentials. 

Remark 3.5. Assumption (A4) is typically verified by solving the Hartree equation in a Sobolev 
space of high index (see e.g. Section l3.2.2[) . Instead of requiring a global-in-time solution tp(-), 
it is enough to have a local-in-time solution on [0, T) for some T > 0. 

Remark 3.6. If sup t 4>(t) < oo, or in other words if and ||y(i)||g 2 are integrable in t over 

R, then all estimates are uniform in time. This describes a scattering regime where the time 
evolution is asymptotically free for large times. Such an integrability condition requires large 
exponents gj, which translates to small exponents pi, i.e. an interaction potential with strong 
decay. 

Remark 3.7. The result easily extends to time-dependent one-particle Hamiltonians h = h(t). 
Replace (Al) and (A2) with 

(Al') The Hamiltonian h(t) is self-adjoint and bounded from below. We assume that there is 
an operator ho ^ that such that ^ h(t) ^ ho for all t. Define the Hilbert space 
Xn = Q(Ei(ho)i) as in (Al). 

(A2') The Hamiltonian Hj^(t) is self-adjoint and bounded from below. We assume that Q{H^{t)) C 
Xn for all t. We also assume that the A r -body propagator U]y(t,s), defined by 

id t U N {t, s) = H N (t)U N {t, s) , U N (s, s) = 1 , 

exists and satisfies E/jv(<, 0)^jv,o £ Q(Hi^(t)) for all t. 
It is then straightforward that Theorem 13.11 holds with the same proof. 

Remark 3.8. In some cases (see e.g. Section [3.2.1 1 below) it is convenient to modify the assump- 
tions as follows. Replace (A3) and (A4) with 

(A3') The interaction potential w is a real and even function satisfying 

||w 2 *M 2 |L < #IM&i (3-6) 
for some constant K > 0. Without loss of generality we assume that K ^ 1. 
(A4') The solution <p(-) of (JOJ) satisfies 
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Then Theorem 13.11 and Corollary 13.21 hold with 



</>{t) = 32K / ds \\<p(s)\\ 



o 



The proof remains virtually unchanged. One replaces (j3.24j) with (|3.6p . as well as (j3.20j) with 

lh*M 2 |L < 2K Mk' 

which is an easy consequence of (13. 6p . 



3.2. Examples. We list two examples of systems satisfying the assumptions of Theorem 13.11 

3.2.1. Particles in a trap. Consider nonrelativistic particles in R 3 confined by a strong trapping 
potential. The particles interact by means of the Coulomb potential: w{x) = A|x| _1 , where 
A £ R. The one-particle Hamiltonian is of the form h = —A+v, where v is a measurable function 
on R 3 . Decompose v into its positive and negative parts: v = u+ — u_, where v+,V- Si 0. We 
assume that v+ £ L\ oc and that V- is — A-form bounded with relative bound less than one, i.e. 
there are constants ^ a < 1 and ^ b < oo such that 

(<p,v-tp) < a(ip, -Atp) + b{ip,tp) . (3.7) 

Thus h + 61 is positive, and it is not hard to see that h is essentially self-adjoint on C^°(R 3 ). 
This follows by density and a standard argument using Riesz's representation theorem to show 
that the equation (h + (b + 1)1)<£ = / has a unique solution ip £ {ip £ L 2 : /i</3 £ L 2 } for each 

/a 2 

It is now easy to see that Assumptions (Al) and (A2) hold with the one-particle Hamiltonian 
h + cl for some c > 0. Let us assume without loss of generality that c = 0. Next, we verify 
Assumptions (A3') and (A4') (see Remark I3.8p . We find 



|2" A 



1^ * M lloo = SU P 



d y i — 12^^)1' 

\x — y\ z 



< (<p,-A<p) < (<p,h<p) + (<p,<p) = \\f\\ 2 Xl 



where the second step follows from Hardy's inequality and translation invariance of A, and the 
third step is a simple consequence of (|3.7|) . This proves (A3'). 

Next, take ipo £ X\. By standard methods (see e.g. the presentation of [7]) one finds that 
(A4') holds. Moreover, the mass ||y(t)|| 2 and the energy 



(<p,h<p) + - I dxdy w(x - y)\(p(x)\' z \(p(y)\'' 



2 



t 

are conserved under time evolution. Using the identity ^ l{|x|^e} e l 2; l~ 2 + l{|x|>e} £ ~ 1 an( i 
Hardy's inequality one sees that 

hmh £ E*{t) + Mt)\\\ 

and therefore ^ C for all t. We conclude: Theorem 13.11 holds with (p(t) = Ct. More 

generally, the preceding discussion holds for interaction potentials w £ L 3 + L°°, where L v w 
denotes the weak L p space (see e.g. [H]). This follows from a short computation using symmetric- 
decreasing rearrangements; we omit further details. This example generalizes the results of |3j, 
P2] and @]. 
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3.2.2. A boson star. Consider semirelativistic particles in M 3 whose one-particle Hamiltonian is 
given by h = Vl — A. The particles interact by means of a Coulomb potential: w(x) = A|x| _1 . 
We impose the condition A > — 4/-7T. This condition is necessary for both the stability of the 
iV-body problem (i.e. Assumption (A2)) and the global well-posedness of the Hartree equation. 
See [Hl[7] for details. It is well known that Assumptions (Al) and (A2) hold in this case. 

In order to show (A4) we need some regularity of <£>(■)■ To this end, let s > 1 and take 
ipo G H s . Theorem 3 of [7] implies that (|1.3|) has a unique global solution in H s . Therefore 
Sobolev's inequality implies that (A4) holds with 

1 1 s 
Vl ~ 2 _ 3' 

Thus qi > 6, and (A3) holds with appropriately chosen values of pi,p2- We conclude: Theorem 
13.11 holds for some continuous function cp(t). (In fact, as shown in [7], one has the bound 
4>{t) < e ct .) This example generalizes the result of pp. 

3.3. Proof of Theorem [Sj} 

3.3.1. A family of projectors. Define the time-dependent projectors 

p(t) := \<p(t))(<p(t)\, q(t) := 1 - p(t) . 

Write 

1 = (Pi + gi) • • • (PN + Qn) (3.8) 

and define Pk, for k = 0, . . . , N, as the term obtained by multiplying out (|3.8p and selecting all 
summands containing k factors q. In other words, 

N 

aGlO.l}^: »=1 

J2i a i= k 

If k 7^ {0, . . . , A^} we set = 0. It is easy to see that the following properties hold: 

(i) Pk is an orthogonal projector, 

(ii) P k p = 5 kl Pk, 

(hi) J2 k p k = i. 

Next, for any function / : {0, . . . , N} — > C we define the operator 

/ := £/(*)JV (3.10) 

k 

It follows immediately that 

fg = fg, 

and that / commutes with pi and P k . We shall often make use of the functions 



m(k) := 
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We have the relation 

^E* = ^EE*^ = ^E^ = ( 3 - n ) 

i k i k 

Thus, by symmetry of ^, we get 

a = = (*,ra*). (3.12) 
The correspondence gi ~ fh of (|3.1ip yields the following useful bounds. 
Lemma 3.9. For any nonnegative function f : {0, . . . ,N} — > [0, oo) we have 

= (*,/ro*), (3.13) 

< -^ T ( l f,/m 2 f}. (3.14) 

Proof. The proof of (|3.13|) is an immediate consequence of (|3.11|) . In order to prove (|3.14|) we 
write, using symmetry of ^ as well as (|3.1ip . 

<*,/9i«2*) = Ar(Ar 1 _ 1) E( lI/ '^g^) 



r (iv-i) 

which is the claim. □ 
Next, we introduce the shift operation r n , n £ Z, defined on functions / through 

(t„/)(*0 := /(fc + n). (3.15) 
Its usefulness for our purposes is encapsulated by the following lemma. 

Lemma 3.10. Let r ^ 1 and ^4 &e an operator on TiyK Let Qi, i = 1, 2, 6e £u>o projectors of the 
form 

Qi = #1 • • • #r > 
where each # stands for either p or q. Then 

Q\A\...rfQl = QlT n fAi_ r Q2, 

where n = n2 — n\ and Hi is the number of factors q in Qi. 
Proof. Define 

N 

p i ■■= e n pi " ■ 

aG{0,l} JV - r i=r+l 
iZi a i= k 

Then, 

Qi? = J2f(k)QiP k = ^f(k)QiP£_ ni = Y,f^ + n i)Qi pr k- 

k k k 

The claim follows from the fact that P£ commutes with A±___ r . □ 



12 



3.3.2. A bound on a. Let us abbreviate 

:= w*\ip\ 2 . 

From (A3) and (A4) we find G L°° (see (13301) below). Then id t ip = (h + W v )ip, where 
h + G >C(Xi; X*). Thus, for any tp G X\ independent of t we have 

id t (4>,p^} = (ip,[h + wv, p ]ip). 

On the other hand, it is easy to see from (A3) and (A4) that m^> G Q(H). Combining these 
observations, and noting that ^ G Q(H) C X by (A2), we see that a is differentiable in t with 
derivative 

a = [H - ,m]*> , 
where IP := ^ + Wf ). Thus, 



d = i( * 



N 

i<j 



By symmetry of \& and m we get 



a = ^(V,[(N-l)Ww-NW?-NW%,rh]V). (3.16) 
In order to estimate the right-hand side, we introduce 

1 = (pi + 9i)(p2 + 92) 

on both sides of the commutator in (|3.16|) . Of the sixteen resulting terms only three different 
types survive: 

i<¥ , PlP2 [(N - 1)W 12 - NWf - NW$ , fh] q lP2 y) (I) 

i<¥ , q lP2 [(N - 1)W 12 - NW? - NW$ , fh] qi q 2 ^>) (II) 

^<tf , PlP2 [(N-l)W 12 - NWf - NW$ ,m] qi q 2 y) (III) . 

Indeed, Lemma 13.101 implies that terms with the same number of factors q on the left and on 
the right vanish. What remains is 

d = 2(1) + 2(11) + (III) + complex conjugate . 

The remainder of the proof consists in estimating each term. 
Term (I). First, we remark that 

P2W 12P2 = P2 Wf . (3.17) 
This is easiest to see using operator kernels (we drop the trivial indices X3, 2/3, . . . , xn, vn)'- 



( P2 Wi 2P2 ){xi,x 2 ;yi,y 2 ) = J dz f(x 2 ) <p(z) w(x\ - z) 5(xi - y\) <p{z) <p(y 2 ) 

= ^p(x 2 )^p(y 2 )S(x 1 - yi) (w * M 2 )(xi) • 
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Therefore, 

(I) = l(y , PlP2 [(N - l)Wf - NWf,m)q lP2 y) = y(*,PiP 2 [Wf ,fh)q lP2 y) . 
Using Lemma 13.101 we find 

(I) = ,pip 2 Wf{m-T^rm)qip2^) = -^{^ ,P\P2Wfqip 2 ^) . 

This gives 

|(i)| < ^ Halloo = ^lh*M 2 IL. 

By (A3), we may write 

w = w^+w®, w {i) E U H . (3.18) 

By Young's inequality, 

\\w {l) * M 2 || < H^H 2 



ioo - „-- wpiWYWn > 
where r\ , r 2 are defined through 

1 = - + -. 3.19 

Pi n 

Therefore, 

lh*M 2 L < Ik^lUlklln+II^IUI^IIr, < (ll^lU + II^Udl^lln + bllr,) 2 - 
Taking the infimum over all decompositions (13.180 yields 

Halloo = ||u;*|^| 2 |L < \\w\\ LPl+LP2 {y\\ ri + \\<p\\ r2 ) 2 . (3.20) 
Note that (A3) and (A4) imply 

2 < n < gi, (3.21) 
so that the right-hand side of (I3.20p is finite. Summarizing, 



1 



|(I)| < 2^11^11^+^ (l^lk + blU) 2 - (3-22) 
Term (II). Applying Lemma 13.101 to (II) yields 

(II) = ^{V,qip 2 {(N-l)W 12 - NW^)(m-r^m)qiq 2 ^) 
= \U , qiP2(^^W 12 -Wn qi q 2 ^Y 

so that 

| (II) | -\(^, qi p 2 Wi 2 qiq 2 ^)\ + -\(^,q lP2 W^qiq 2 ^)\. (3.23) 
The second term of (|3.23p is bounded by 

J||PnUI|gi*|| 2 < h\w\\ LP1+LP2 {M ri + M r2 ) 2 a, 



14 



where we used the bound (|3.20p as well as (|3.12p . 

The first term of (|3.23j) is bounded using Cauchy-Schwarz by 

2 QiP2Wf 2 p 2 qi^)V(^ , <M2*> = l\J(^,qiP2(w 2 * \<p\ 2 ) ^qi^Vi^'Qi Q^) 

This follows by applying (13. 1TH to W 2 . Thus we get the bound 

^Hgi^llVlh 2 * M 2 IL = ^ a \J\\ w2 * M 2 IL- 

We now proceed as above. Using the decomposition (|3.18j) we get 

K*M 2 |L < 2||( U ;W) 2 *M 2 |L + 2||( U ;( 2 )) 2 *M 2 |L- 
Then Young's inequality gives 



||(ru») 2 *M 5 



< \\w 



(0|| 2 11-112 
I Pi 1 



m\ qi 



which implies that 



\\w 2 * IvHL < 2 \H\LPi+LP2(\M qi + IMI<? 2 ) 2 ■ 
Putting all of this together we get 

|( n )| < ^\\w\\lpi+lp2 \/2(|M| ?1 + H^llga) + (iMIn + WvWn)' 



(3.24) 



a . 



Term (III). The final term (III) is equal to 
7;(^,PiP2[(N - l)Wi2,m] qi q 2 y) = ^{y,PiP2(N-l)Wi 2 (m-T^fn)q 1 q 2 i Sf) 



i N N 1 ,PlP2W 12 qiq2^) , 



where we used Lemma 13.101 Next, we note that, on the range of q±, the operator n 1 is well- 
defined and bounded. Thus (III) is equal to 

N - 1 



* ,pip 2 Wi 2 nn qiq 2 ^) 



l ^-, ^ ,pip 2 T^nWi 2 n 1 q\q 2 fy) 



N x L r 3 " ' N 

where we used Lemma 13.101 again. We now use Cauchy-Schwarz to get 



| (HI) | < 



,pip 2 T 2 nW? 2 r 2 npip 2 ^)J ($,fi 2 qiq 2 ^) 



< 



^ ,PiP2TYn (w 2 * \(f\ 2 ) 1 T 2 n p\p 2 ^>) y (\l/,m l qiq 2 ^>) 



w 2 * \(p\ 2 



w 2 * 



w 2 * 



w 2 * 



N 



N - 1 



N 



N - 1 

N 



<p\' 



N 



<P\ : 



N - 1 

N 
N - 1 



I2a\ 



2 [a H 
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Using the estimate (|3.24p we get finally 

| (III) | ^ 2V2\\w\\ LPl+LP2 (y\\ qi + \\ip\\ q2 ) y 



N 



a+ N 



Conclusion of the proof. We have shown that the estimate (|3.2p holds with 



BN(t) = 2\\w\\lpi+LP2 
A N {t) - 



(|b(i)||n + lb(*)l|r 2 ) 2 + 6(||^)|| gi + 



Using L 2 -norm conservation ||<^(t)|| = 1 and interpolation we find 

B N (t) < 16\\w\\ LPl+LP2 (\\<p(t)\\ qi + y(t)\\ q2 ) . 
The claim now follows from the Gronwall estimate (13.31). 
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Thus, 



4. Convergence for stronger singularities 

In this section we extend the results of the Section to more singular interaction potentials. We 
consider the case w £ L Po + L°° , where 

Po 2d 

For example in three dimensions po = 6/5, which corresponds to singularities up to, but not 
including, the type |x| _5//2 . Of course, there are other restrictions on the interaction potential 
which ensure the stability of the iV-body Hamiltonian and the well-posedness of the Hartree 
equation. In practice, it is often these latter restrictions that determine the class of allowed 
singularities. 

In the words of [TT] (p. 169), it is "venerable physical folklore" that an iV-body Hamiltonian 
of the form (|3.4p . with h = —A and w(x) = \x\~^ for £ < 2, produces reasonable quantum 
dynamics in three dimensions. Mathematically, this means that such a Hamiltonian is self- 
adjoint; this is a well-known result (see e.g. [H])- The corresponding Hartree equation is known 
to be globally well-posed (see [5]). This section answers (affirmatively) the question whether, in 
the case of such singular interaction potentials, the mean-field limit of the iV-body dynamics is 
governed by the Hartree equation. 

4.1. Outline and main result. As in Section [3l we need to control expressions of the form 
\\ w2 * M 2 ||oo- The situation is considerably more involved when w 2 is not locally integrable. An 
important step in dealing with such potentials in our proof is to express w as the divergence 
of a vector field f £ L 2 . This approach requires the control of not only a = \\qi^f\\ 2 but also 
llViqi^H 2 , which arises from integrating by parts in expressions containing the factor V • £. As 
it turns out, (3, defined through 

0N(t) := (V N ,n* N )\., (4.2) 
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does the trick. This follows from an estimate exploiting conservation of energy (see Lemma 14.61 
below). The inequality m ^ n and the representation (I3.12|) yield 

a sC 0. (4.3) 

We consider a Hamiltonian of the form (13.4|) and make the following assumptions. 

(Bl) The one-particle Hamiltonian h is self-adjoint and bounded from below. Without loss of 
generality we assume that h 0. We also assume that there are constants K\, k 2 > such 
that 

-A < K\ h + k 2 , 

as an inequality of forms on Tiy-K 

(B2) The Hamiltonian (|3.4p is self-adjoint and bounded from below. We also assume that 
Q(Hn) C Xjy, where Xn is defined as in Assumption (Al). 

(B3) There is a constant K3 G (0, 1) such that 

< (1- K 3 )(h 1 + h 2 ) + W 12 , 

as an inequality of forms on TL^ 2 \ 

(B4) The interaction potential if is a real and even function satisfying w G L p + L°°, where 
po < p ^ 2. 

(B5) The solution (/?(•) of (jl.3j) satisfies 

p(-) G C(M;X 1 2 nL 0O )nC 1 (M;L 2 ), 
where X 2 := Q(h 2 ) C L 2 is equipped with the norm 

|M| X? := \\{l + h 2 f' 2 tp\\. 

Next, we define the microscopic energy per particle 

E%(t) := ±(* N ,H N * N )\ t , 

as well as the Hartree energy 

E?{t) := {<p,h<p) + ± J dxdy^(x-y)|^(x)| 2 |^(y)| 2 

By spectral calculus, E^(t) is independent oft. Also, invoking Assumption (B5) to differentiate 
E^(t) with respect to t shows that E^{t) is conserved as well. Summarizing, 

E%{t) = E%(0), E?(t) = EV(0), teR. 

We may now state the main result of this section. 
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Theorem 4.1. Let ^n,o G Q{Hn) o,nd assume that Assumptions (Bl) - (B5) hold. Then there 
is a constant K , depending only on d, h, w and p, such that 



Mt) < (p N (0)+E N -E^ + ^je K ^, 



where 



and 



p/po - 1 



2p/p -p/2-1 



(4.4) 



4>{t) := / ds 1 + IMOH 



\x%nL°° 



Remark 4.2. We have convergence to the mean-field limit whenever limjvi?* = E v and 
limjv /3jv(0) = 0. For instance if we start in a fully factorized state, ^jv,o = <Po N i then (3m (0) = 
and 



E% - E^ = — ((f (g> ip , W12 <fo ® <po) 



so that the Theorem 14. 1 1 yields 

and the analogue of Corollary 13.21 holds. 

Remark 4.3. The following graph shows the dependence of rj on p for d = 3, i.e. po = 6/5. 




P 



Remark 4.4. Theorem 14.11 remains valid for a large class of time-dependent one-particle Hamil- 
tonians h(t). See Section [4.41 below for a full discussion. 



Remark 4.5. In three dimensions Assumption (Bl) and Sobolev's inequality imply that \\(p\ 
\\<p\\xf, so that Assumption (B5) is equivalent to ip G C(R;Xl) n C^MjL 2 ). 



< 
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4.2. Example: nonrelativistic particles with interaction potential of critical type. Consider 
nonrelativistic particles in R 3 with one-particle Hamiltonian h = —A. The interaction potential 
is given by w(x) = A|x| -2 . This corresponds to a critical nonlinearity of the Hartree equation. 
We require that A > —1/2, which ensures that the A^-body Hamiltonian is stable and the Hartree 
equation has global solutions. To see this, recall Hardy's inequality in three dimensions, 

{(p,\x\~ 2 <p) ^ 4{<p,-A<p). (4.5) 

One easily infers that Assumptions (Bl) - (B3) hold. Moreover, Assumption (B4) holds for any 
p < 3/2. 

In order to verify Assumption (B5) we refer to [5], where local well-posedness is proven. 
Global existence follows by standard methods using conservation of the mass \\f\\ 2 , conservation 
of the energy E^, and Hardy's inequality (|4.5p . Together they yield an a-priori bound on ||y||xi) 
from which an a-priori bound for ||y||x 2 may be inferred; see [5] for details. 

We conclude: For any rj < 1/3 there is a continuous function <j>(t) such that Theorem 14.11 
holds. 

4.3. Proof of Theorem 14.11 

4.3.1. An energy estimate. In the first step of our proof we exploit conservation of energy to 
derive an estimate on || Vi^i^H. 

Lemma 4.6. Assume that Assumptions (Bl) - (B5) hold. Then 

HVigi^H 2 < E*-E*+(l + \\<p\\ 2 x , nLoo )(p + 

Proof. Write 



as well as 



Inserting 



E? = (<p,htp) + ±(<p,W*<p), (4.6) 
E* = (y,h 1 V) + ] -^{y,W 12 y). (4.7) 



1 = PlP2 + (1 -PlP2) 

in front of every in (|4.7p and multiplying everything out yields 

(¥,(l-piJ>2)Ml-PlP2)tf) 

= E* — (^,PlP2hiPip 2 ^) 

N - 1 

,PlP2W 12 piP2^} 

- (^,(l- Pl p2)h lPl p2^) - (*,PU>2/»l(l-piP2)*> 

" -^^(^,(1-P1P2)W12P1P2^) " ^^(^ ,PlP2W U {l - pm)^) 
~ ^^(^,{^-PlP2)Wi 2 (l-piP2)^) ■ 
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We want to find an upper bound for the left-hand side. In order to control the last term on the 
right-hand side for negative interaction potentials, we need to use some of the kinetic energy on 
the left-hand side. To this end, we split the left-hand side by multiplying it with 1 = K3 + (l — K3). 
Thus, using (14. 6ft . we get 



k 3 (¥,(1 -Pip 2 )hx(l -pip 2 )V) 

= E^-E^ 

- (* ,PiP2hipiP2^) + {p>,hip) 
N 



^(^, PlP2 W 12 PlP2^) + -(ip^if) 



(^,(l- PlP 2)h lPl p2^) 

N - 1 



2N 
N - 1 

2N 



* , (1 - PIP2)W\2P1P2^ 



N — 1 



2N 



(tf ,PlP2Wi 2 (l -P1P2)*) 



-pip 2 )W 12 (l ~PiP2)^) - (1 - Ks)<*,(l-PiP2)^i(l -PiPa)*> • (4- 



The rest of the proof consists in estimating each line on the right-hand side of (|4.8p separately. 
There is nothing to be done with the first line. 

Line 6. The last line of (j4.8jl is equal to 



-gjy^ (* » (! - P1P2)W 12 (1 ~ PU>2)*) - ^(1 - K 3 )(* , (1 - PlP2)0l + /l 2 )(l " P1P2)*) 



JV- 1 
2iV 



* , (l-piP2)[(l-K3)(^l + ^2)+Wi2](l-piP2)^) < 



where in the last step we used Assumption (B3). 
Line 2. The second line on the right-hand side of 



is bounded in absolute value by 



\((f,h(p) - ,piP2hiPiP2^)\ 



= ((p,fup)\(V,(l -pi|>2)*)| 

= (p, hip) |(#, (qxp2 + pi q 2 + gi^)*) 

^ 3a (p, hip) 

^ 3/3(ip,hip), 



where in the last step we used (|4.3p . 

Line 5. The third line on the right-hand side of 



1 2V- 
-b.P^w 



(*,PlP2Wi 2 PlP2*) 



1 



is bounded in absolute value by 
N- 1 



1 



N 



■(*,PlP2*) 



< 2 ll ^ l|c 



< 2 ll ^ l|c 



1 



(*, (gip 2 + Pl<72 + gi^)*) + J^(^ ,PlP2^) 



1 

TV 
1 

iV 



a + 
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As in (|3.20p . one finds that 

Halloo < IMIli+L-IM^^oc • 

Line 4- The fourth line on the right-hand side of (|4.8p is bounded in absolute value by 

|(#, (1 - piP2)hipip2^)\ = \(^,(qiP2+piq2 + qiq2)hipiP2^)\ 

= \(^,qih lP ip2^)\ 

= ,qin~ 1/2 n 1/2 hipip 2 $>)\ 

I i 1 /q ... -— 1/2 \ I 

= \{fy,qin '^hiTin 1 p\P2^)\ , 
where in the last step we used Lemma 13.101 Using Cauchy-Schwarz, we thus get 



- PiP2)h 1 p 1 p 2 ^)\ < J^^q^ 1 ^) J(^, Pl p 2 mi 1/2 h 2 ^n /2 Pl p 2 ^ 



= ^{^,n^)^{i P ,h 2 ip)^(^>,Trfip 1 p2^), 

where in the second step we used Lemma [3791 Using 



(rm)(fc) = J^±I < n(k) + 



we find 



|<*,(1 - pip 2 )hipi P 2^) | s£ ^fp^/(tp,h 2 ^)^(^ ,n^>) + ^ 

^ 2y/(<p,h* V )(p+-j=\ 



. N 

ATl/4 



Line 5. Finally, we turn our attention to the fifth line on the right-hand side of (|4.8p . which is 
bounded in absolute value by 

\(*,PiP2W 12 (l -P1P2)*>| = \(^,PiP2Wi 2 (piq2 + qiP2 + qiq2^>)\ < 2(a) + (b) , 

where 

(a) := \(^ ,p\P2Wi2q\P2^)\ , (b) := | piP2Wi 2 qi qi®) \ ■ 
One finds, using (|3.17j) . Lemma f3. 101 and Lemma 13791 

(a) = |(*, Pl p 2 Wfgi*)| 

= |<* ,pip 2 Wfn 1/2 n~ l/2 qi®)\ 
= |<* ,pip 2 nn ,2 Wfn~ 112 gi*)| 



^ 2||W|| 00 f/? + 
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The estimation of (b) requires a little more effort. We start by splitting 



w 



+ w 



(oo) 



w 



^ G LP , G L c 



This yields (b) ^ (b)^ + (b)^ 00,1 in self-explanatory notation. Let us first concentrate on (b) 



(oo). 



(b) 



(oo) 



{m, Pl p 2 wlf qi q 2 y)\ 

,p\p 2 w[ 2 nrT 1 qiq 2 ^)\ 
, pip 2 T2~n wffl n~ l qiq 2 ^>) 



l^ (oo) HooV<^,^ 2 ^> J(V,ft- 2 qiq 2 y 



w oo 



a + N 



< 21 to 



Mi 



N 



Let us now consider (b)^ p \ In order to deal with the singularities in we write it as the 
divergence of a vector field £, 



(p) 



(4.9) 



This is nothing but a problem of electrostatics, which is solved by 



c JUL * 



with some constant C depending on d. By the Hardy-Littlewood-Sobolev inequality, we find 

1 1 1 



w 



(p)\ 



(4.10) 



p q p d 

Thus if p ^ po then q ^ 2. Denote by X12 multiplication by — x 2 ). For the following it is 
convenient to write V • £ = V p £ p , where a summation over p = 1, . . . , e£ is implied. 
Recalling Lemma 13.101 we therefore get 



(b) 



(p) 



= |<*,mP2T 2 n(V?X' , )i2n 1 q 1 q 2 V)\. 

Integrating by parts yields 

(b) (p) ^ \{V{vmpip 2 ^,Xl 2 rr 1 q 1 q 2 V)\ + \{^ Pl p 2 ^ ,X^ 2 V[n- l qiq 2 ^)\ 



(4.11) 



Let us begin by estimating the first term. Recalling that p = \ip)(ip\, we find that the first term 
on the right-hand side of (|4.1ip is equal to 



\{X^ 2 p 2 {V p p)iT 2 n^ ,ri l qiq 2 $?) 



< 



(VPp)i^^ ,p 2 X p 2 X^ 2 p 2 (V a p)iT 2 ~n^) \\n V<?2^| 



< \/llM 2 *£ 2 IL llv^H 11^^11 \\n 1 qmM 



^ U\U \\<p\\L*nL°° WfWxiX/oi + — Va, 
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where we used Young's inequality, Assumption (Bl), and Lemma 13.91 Recalling that [3 ^ a, we 
conclude that the first term on the right-hand side of (14. lip is bounded by 

cy\\ 2 Xl nL~ (P+Jj)- 

Next, we estimate the second term on the right-hand side of (14.111) . It is equal to 



[X^ 2 pip2T2n^ ,Vln l qiq2^>)\ ^ \J {T2n * , PiP2X'( 2 pip2 T2~n *) || Vi n 1 qiq2^\ 



2 *C 2 |L ll^n^ll ||Vin 1 gig2*| 



< ll£lUMlL 2 nL°°y « + ^ ||Vin 1 gig 2 *||. 

We estimate || Vi n _1 91 g^*)! by introducing 1 = p\ + q\ on the left. The term arising from p\ 
is bounded by 

llpiVin" 1 qi 92*11 = \\p1q2T1n~ 1 Vi9i*|| 



/ 1 - 

\ \ Viqi ^ ' jrriJ2i iTin 2 Vi9i* 



N 



1 - 

V191* , — ^qiT^H 2 Vi9i^ 



i=i 



n 2 rTn 2 V191*) 
< ||Vi9i*||. 

The term arising from 91 in the above splitting is dealt with in exactly the same way. Thus we 
have proven that the second term on the right-hand side of (|4.1ip is bounded by 



/5 + ^||Vi9i*| 



Summarizing, we have 



(b) (p) < IMft in L~ (P+ ^) + \W\\l^l~\JJ+^ l|Vi9i*| 

Conclusion of the proof. Putting all the estimates of the right-hand side of 
find 

(*,(1-PU>2)^1(1-P1P2)*> 



together, we 



■ H*?nL~) (P+ 7=) + IMUanL-y/9+ ^ l|Vi9i*|| . (4.12) 
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Next, from 1 — P1P2 = Piqi + qi we deduce 

llV^i^H = ||v^i(l-PiP2)*- V^iPi92*|| < ||\Ai(l-PiP2)*|| + ||V^iPig2*||- 
Now, recalling that p = \<p)(<p\, we find 

||-nAipi«2*|| < H-nAipiIII|92*|| < WvWxiVP- 

Therefore, 



^i9i*H 2 £ ||\Ai(l-piP2)*|| 2 + IMIx^ • 
Plugging in (|4.12|> yields 

||V^?i*|| 2 < - ^ + (1 + Ikll^nicc) (/? + -!=) + \\<p\\u>nL<*>)Jp + ^ l|Vi?i*|| . 
Next, we observe that Assumption (Bl) implies 

||Vigi*|| < ||v^i?i*|| + y/P, 

so that we get 

ll^i^H 2 < - E* + (1 + |M& ?nLO o) (/? + + IMUw^ + ^ \\Vh~mn ■ 

Now we claim that 

iiv^«i*ir s E^-E^+(i+y\\^ nLX )^+-^=y (4.i3) 

This follows from the general estimate 

x 2 < C^ + ax) =► x 2 < 20R + CV, 
which itself follows from the elementary inequality 

C(R + ax) < CR+^C 2 a 2 + ^x 2 . 

The claim of the Lemma now follows from (|4.13p by using Assumption (Bl). □ 

4.3.2. A bound on j3. We start exactly as in Section [3l Assumptions (Bl) - (B5) imply that (3 
is differentiable in t with derivative 

p = ^(y,[(N-l)W 12 -NW¥ -NW£,n]V) 

= 2(1) + 2(11) + (III) + complex conjugate , (4.14) 
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where 

(I) := ^(y,p lP2 [(N-l)W 12 -NW? -NWg,n]q lP2 V) 

(II) := ,qi P2 [(N -l)W 12 - NWf - NW$ ,n] qi q 2 ^) 

(III) := i(*,p lP2 [(jV-l)W 12 - NWf- NW$,n]qtf2*) 

Term (I). Using (|3.17p we find 

2 1 (I) | = \{V,pip 2 [(N -1)W 12 - NWf - NW? ,n]q lP2 V) 
= \{y,PiP2[W¥,n]qi P2 y)\ 
= \{&,PiP2W?(n- T^n)qxp 2 ^)\ , 

where we used Lemma 13.101 Define 



M (fc) := N(n(k) - (r_m)(fc)) = ^ n^) , fc = 1, 



Thus, 



| (I) | = ±\(9,piP2W?jlqiP2*)\ 
1 



~ j^llyllianLco , 

by (|513j> . 

Term (II). Using Lemma 13.101 we find 

2|(II)| = |<tt,gij*[(tf-l)Wi 2 - | 
= , 9lP2 f ^^Wia - W-f J J2qiq 2 ^ 

^ 1<* , gip 2 TUi 2 ju gig 2 * ) | + | , giP2Wf A» 9192*) 



= : (a) =: (b) 

One immediately finds 

(b) < n^iiooiki^iiv^,^!^} < y\\l* nL ~p. 

In (a) we split 

with a resulting splitting (a) ^ (a) + ( a ) ■ The easy part is 

(a)(°°) < H^IUHgx*!! 2 < /?. 
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In order to deal with (a)^ we write = V • £ as the divergence of a vector field £, exactly as 
in the proof of Lemma l4,6t see (14. 9\\ and the remarks after it. We integrate by parts to find 



(a) 



(P) = |<*, 9 iP2(V?X")l2At9l92*)| 



(4.19) 



The first term of (I4.19P is equal to 



\{X p 12 p 2 V{ qi ^ ^ qi q 2 ^)\ ^ J(V p 1 q 1 *,P2Xe 2 X? 2 p 2 V° 1 q 1 *)J(*,Fq 1 q 2 * 



< ^lle 2 *^| 2 ||oo||V igi M/|| 



< ^ll,lkllLanL»||Vi?i*||v^ 

< l|Vigi^|| 2 IM^m/- + P llylU 2 ni- . 



where in the second step we used (14.15p . in the third Lemma 13.91 and in the last (14, 3h . Young's 
inequality, and (|4,10p . The second term of (|4.19p is equal to 

| (q lP2 V , X P 2 ( Pl + qi )V{ ju qi q 2 ^>) | 

< |(giP2*, X P 2Pl TnLV p qi q 2 y) | + \(q lP2 ^,X P 2qi ^V?gig 2 *)| , (4.20) 

where we used Lemma [3. 101 We estimate the first term of (I4.20p . The second term is dealt with 
in exactly the same way. We find 



( Pl X P 2qiP2 ^ ,T!fMV^qiq 2 ^)\ < J(y,qi P2 Xf 2P2 q 1 y)J(V 1 q 1 y,q 2 nfiq 2 V 1 q 1 y 



< VU 2 * M 2 Hoo y(Viqi^,n~ 2 q 2 Viqi^) 

~1 * 

HCllg IMU 2 nZ°° V 7 " a N _ 7 Yl < V i9i^ , n~ 2 q&iqi 

\ i=2 



1 * 

^ i=i 



iV 
iV- 1 



Vigi^, n~ 2 n 2 ViQi*) 



< \M 



MlL 2 nL°°V P 

v^IIVi?i*|| 

< /3 IMIjy*n£°° + llVigi^H 2 |M|i2 nL oo . 

In summary, we have proven that 

|( n )| £ 0Hl£3n£<» + ||Vigi*|| 2 H^H^nioo . 

Term (III). Using Lemma 13.101 we find 

2|(IH)| = (AT- l)|(*,pip2 [W 12 ,n] gi q 2 V)\ = (N - , P ip 2 W 12 (n - f^n) qi q 2 ^) 
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Defining 

u(k) := N(n(k)-(T_ 2 n)(k)) = - ^L— < n'^k) , k = 2,...,N, (4.21) 

\l k + \J k — I 

we have 

2|(III)| s£ \{^ ,PiP2W 12 v qiq2^)\ 

As usual we start by splitting 

w = W (P) + w i°°) , W <P) e LP , G L°° , 

with the induced splitting (III) = (III) (p) + (III) (oo) . Thus, using Lemma 13.101 we find 
2|(III) (oo) | = | , vwMf nV2 ^-V2 p qiq2 ^ | 

= | (* > PlP2 ^l ( 2° 0) ^" V2 ? 9l92* ) | 

,Hl 



where in the fifth step we used Lemma 13.91 

In order to estimate (III)^ we introduce a splitting of into "singular" and "regular" 
parts, 

w (p) = W M + >,2) := i {|wW|>a} + >) l {|wW|<a} , (4.22) 

where a is a positive (iV-dependent) constant we choose later. For future reference we record 
the estimates 

lM P,1) IU < a 1-p/P0 Ik^H^ , (4.23a) 
||>' 2) || 2 ^ a^Ww^Wf. (4.23b) 

The proof of (|4.23p is elementary; for instance (|4.23ap follows from 
\\ w M\\Po = [dx \wV>\ p \wto\ p °- p l, 



^ a P0 ~ p I dx |w (p) | P l { | tt (p)| >a} ^ a P0 " p / dx \w^\ p . 



Let us start with (III)^' 1 ^. As in (|4.9p . we use the representation 

w ( P ,i) = v-£. 
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Then (EL"TOl) and (I4.23al) imply that 

Uh < \\™ M \\ Po < ^ Plm - (4-24) 

Integrating by parts, we find 

2|(III)^ 1 )| = \(*, P i P 2W& 1) vq 1 q 2 *)\ 

= \(^, P m(V p 1 X^)uq 1 q 2 ^)\ 

< \{V p lPl p 2 ^,X p 2 uq 1 q 2 yS)\ + \(p 1 p 2 ^,X P 2 V p 1 uq 1 q 2 ^)\. (4.25) 
Using ||Vp|| = ||Vy|| and Lemma 13, 91 we find that the first term of (|4.25j) is bounded by 

^(v?pi* ,P2A? 2 xf 2 P2v?pi*> , ^ qiq2 ^) < iivpini^iuiieibv^ 

< HV^IIblloca 1 -^ v 7 ^ 

< llV^IHI^IU^ + a 2 - 2 ^), 

where in the second step we used the estimate (14.24j) . Next, using Lemma [3.10| we find that the 
second term of ()4.25p is equal to 

\(p lP 2^,X P 2 (p 1 +q 1 )V p Dq 1 q 2 ^)\ 

sC \{pip 2 ^,X P 2Pl f^V^qiq 2 ^)\ + | (pip 2 ^, X^ qi vV{qi q 2 ^)\ ■ 

We estimate the first term (the second is dealt with in exactly the same way): 

\(p lP2 ^ ,X p 2Pl n^V p qi q 2 ^)\ sC yj (^, Pl p 2 Xf 2Pl p 2 ^) yj (Vift^ , n^ 2 q 2 Viqi^) 



I 1 

< \J\\P2 X 12P2L j^—^J2(Viqiy,n~ 2 qiViqi* 

\ i=2 

\ i=i 



^ II^IU^-^ + IIVx^n 2 ). 

Summarizing, 

< I^IU^II^ + IIVigi^f + a 2 - 2 ^ !!^). 

Finally, we estimate 

(III) (p ' 2) = \(^,pip 2 W^ 2) uq iq2 ^)\ = \(y, PlP 2Wif +^)) qi q 2 ^)\, (4.26) 



where 



1 = X W +X i2 \ X (1) ,X (2) € {0,1}{°<"-^ 
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is some partition of the unity to be chosen later. The need for this partitioning will soon become 
clear. In order to bound the term with we note that the operator norm of pip 2 W^ q\q2 
on the full space L 2 (M, dN ) is much larger than on its symmetric subspace. Thus, as a first step, 



we symmetrize the operator Q1Q2 i n coordinate 2. We get the bound 

\{^,p lP2 W^ 2) u^m qi q 2 ^)\ 



1 



N — 1 
1 

N-l 



N 



i=2 



* . J2pmwlf' 2) mi x (1) £gi* 
i 

vqx$ 



A' 



]T (V^lPiW^qtqixWqiqMj %P&) 

i,j=2 



Using 



we find 



where 



||z?<7i'I>|| ^ \\n 1 qi i &\\ < 1 



A 



£ , PlPiW^ 2} q iqi X « qMjPjP^) 

N 

B := {^,PlPiW^ 2) q iqi X« w[f ] p lPl ^) . 



i=2 



The easy part is 



A' 



5 < Y,(y,piPi(wfr 2) ) 2 PiPiy 



i=2 
N 



II 

< (^-i)ll^llLll>' 2) lll 



< Na 2 -P\\ V \ 



2 

00 



Let us therefore concentrate on 



= + A 2 , 



(4.27) 
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with A = A\ + A2 arising from the splitting qi = 1 — p\. We start with 



I Ai| < ]T \(*,PiPi<lj : ^wfr 2) W$ 2): ^q iPj piy)\ 



by Cauchy-Schwarz and symmetry of VP. Here \/- is any complex square root. 
In order to estimate this we claim that, for i 7^ j, 



PiPi\W^ 2) \\W^ 2) \p lPi 



<» 2 >|*M 2 ||* 



(4.28) 



Indeed, by (|3.17p . we have 

The operator pi (|u/ p ' 2 ) | * |(^| 2 ) 1 1 W^' 2 ^ \p± is equal to fjPx, where 

f( Xj ) = [ dx 1 ^^{\w^\*\^\ 2 )(x 1 )\w^(x 1 -x j )\i P (x 1 ). 



Thus, 



Ml * l,ol 2 l| 2 
1 00 



< ||K P ' 2) | * \<p\ : 



from which (|4.28p follows immediately. 
Using (|4.28p . we get 



^ (p ' 2) l*l^ 2 ||L|k 2 x (1) gi^|| 2 



< iV 2 ||>)|| 2 ||^ 2n ^<^,SWg^ 

< iY 2 ||< 2nL ^<*,So)n 2 *). 



Now let us choose 



for some 5 £ (0, 1). Then 
implies 



X«(fc) := l {fc<iV i-n 
(r 2X (1) )n 2 < AT 5 



l^il < [I ^[lianioo JV 



2-5 



(4.29) 
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Similarly, we find 

< 2 H>> 2 )*iH 2 llL(*'S«9i*} 



< N 



2||,„||4 



L 2 nL°° 



L 2 nL°° 



iV 



2-5 



Thus we have proven 



14 < IMlV-^ 2 " 5 



Going back to (|4.27|) . we see that 

|<*,PiP2Wg' 2) ^JW (ZKfe*) | < Hvll^niooJV-*/ 2 + HvllooJV-Vaa 1 -^ 2 
What remains is to estimate is the term of (III)^' 2 ^ containing 

N 



1 



N - 1 



N - 1 



i=2 



N 



^2 (V,PiPiW$®qiqixWvqiqjW}!} 2) PjPi^) ■ 

\ i,3=2 



Using 
we find 

where 



?V2 



2<i^'<AT 
AT 

B := ^<*,PiftW 1 ^ 2) ? i%x (2) ^W 1 ( f' 2 Wi*>- 



(4.30) 



i=2 



Since 
we find 



X {2) {k) = l {k>N i-s } 



X {2) v < X (2) n- X < iV 5 / 2 . 
Thus, ||gi%x (2) ?ll < N 5 / 2 and we get 

b < iy^^^^ip^^) 2 ^*) < iv 1+5 / 2 ||(>- 2 )) 2 *M 2 | 



i=2 
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by (I423E1) . 

Next, using Lemma l3,10( we find 

= £ ^^ 1/2 ^ ,3) 91^ ,3) ^^ 1/2 W1*> 

= ^1+^2, 

where, as above, the splitting A = A± + A 2 arises from writing gi = 1 — p±. Thus, 

\A,\ < £ K^Pl^^T^ 1 / 2 ^'^^^^^^ 
2<^j<iV 

by Cauchy-Schwarz and symmetry of \&. Using (|4.28p we get 

< N 2 Mi 2nLOC (3. 

Similarly, 



\A 2 \ < E I <^ > P^i Pi^i ( f ' 2 Vi^ ( f 2) pi ^ 1/2 I 

< iV 2 ||>)|| 2 |K 2nLoo (M/,n*) 

< iV 2 |K 2nLOO /3. 

Plugging all this back into (|4.30p . we find that 

\(*,p 1 p 2 wi p 2 > 2) v^ qiq2 n>)\ < /?(|M|£ 2nLOO + Halloo) + ll^Hooa 2 -^ 2 - 1 

Summarizing: 

|(III)^| < (l+|| (/ ,||2 2nLoo )^ + a 2- PiV <5/2-l + N S/2 +N -l/2 a l-p/2j t 
from which we deduce 

|(III) (P) | < ll^lloollV^^H 2 

+ (1 + \\<p\\ XinLao ) f/3 + a 2 ~ p N s ' 2 ~ l + N~ s ' 2 + N'^a 1 ^ 2 + a 2 ~ W 
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Let us set a = ajy = i\K and optimize in 8 and £. This yields the relations 



C0-,) + « = i. -5 = «c(i-^ 

which imply 

5 y/po - 1 



2 2p/p - p/2 - 1 ' 

with (5^1. Thus, 

< |M| 00 ||Vi< ? i*|| 2 + (l + ||^||x 1 nL-)(/3 + iV^ 



where rj = 8/2 satisfies 
Conclusion of the proof. We have shown that 

£ ll^llL2nL-l|Vigi^|| 2 +(l + ||^||x 1 nL-)(/3 + iV" ,? ). 
Using Lemma 14.61 we find 

< (l + IMI^nL-) (/* + ^ - & + ^) ■ (4-31) 
The claim then follows from the Gronwall estimate (13.31). 



4.4. A remark on time-dependent external potentials. Theorem 14.11 can be extended to time- 
dependent external potentials h(t) without too much sweat. The only complication is that 
energy is no longer conserved. We overcome this problem by observing that, while the energies 
E^ it) and E^(t) exhibit large variations in t, their difference remains small. In the following 
we estimate the quantity E (t) — E v it) by controlling its time derivative. 

We need the following assumptions, which replace Assumptions (Bl) - (B3). 



(Bl') The Hamiltonian hit) is self-adjoint and bounded from below. We assume that there is 
an operator ho ^ that such that ^ h(t) ^ ho for all t. We define the Hilbert space 
Xn = QwZi{ho)i) as in (Al), and the space X\ = Q(/iq) as i n (B5) using ho- We also 
assume that there are time-independent constants ki,K2 > such that 

-A < K\ hit) + K 2 

for all t. 

We make the following assumptions on the differentiability of h(t). The map t >— > {ip , h{t)ip) 
is continuously differentiable for all ip G Xi, with derivative {tp , h(t)ip) for some self-adjoint 
operator h(t). Moreover, we assume that the quantities 

<¥>(*) , h{t)\{t)) , ||(1 + hit))" 1 / 2 hit) (1 + Ht))- 1 ^ 

are continuous and finite for all t. 
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(B2') The Hamiltonian Hj^(t) is self-adjoint and bounded from below. We assume that Q(H^(t)) C 
Xn for all t. We also assume that the iV-body propagator U]y(t,s), defined by 

id t U N (t, s) = H N (t)U N (t, s) , U N (s, s) = t, 

exists and satisfies U^(t, O^a^o 6 Q{HN(t)) for all t. 

(B3') There is a time-independent constant E (0, 1) such that 

< {l-^){hi{t) + h 2 {t)) + W l2 

for all t. 

Theorem 4.7. Assume that Assumptions (Bl'j - (B3'), (B4), and (B5) hold. Then there is a 
continuous nonnegative function eft, independent of N and *S>n,o> such that 



with 7] defined in (|4.4p . 

Proof. We start by deriving an upper bound on the energy difference £(t) := £7* (t) — E^it). 
Assumptions (Bl') and (B2') and the fundamental theorem of calculus imply 

8{t) = 8(0) + f ds ( {*(s) , kjs^js)) - (<p(s) , h(s)tp(s))) ■ 

= :G(s) 

By inserting 1 = pi(s) + qi(s) on both sides of hi(s) we get (omitting the time argument s) 

G = {W,p 1 h 1 px*)-{<p,h(p) + 2Re{*,p 1 h 1 q 1 V) + {V,q 1 h 1 q x V). (4.32) 
The first two terms of (I4.32j) are equal to 

((*,Pi*> -l)((p,h<p) = a(<p,h(p) < p\((p,htp)\. 
The third term of (|4.32p is bounded, using Lemmas 13.91 and 13.101 by 

2|<^,p 1 A 1 n 1/2 n- 1/2 g 1 ^}| = 2\ (h lPl rrh 1/2 $ , n~ 1/2 qi$) \ 



< y(^ 1/2 *,pi/ifpirTn 1/2 ^) ||ri _1/2 gi^|| 

< yj\{<p,h 2 <p)\ ^(^,nn~y) y/^.w- 1 ?!^) 

< J\(< p ,h*<p)\(() + -L 
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The last term of f|4.32j) is equal to 



+ hi) 1/2 (l + hy^h^l + /ii)- 1/2 (l + hf^q^) 

< ||(1 + /i)" 1 / 2 ^! + /i)" 1/2 || ||(1 + /ii) 1/2 gi^|| 2 . 

Thus, using Assumption (Bl') we conclude that 

G(t) < C{t)Ut) + -L + ||M*) V V(W)f) (4-33) 



for all t. Here, and in the following, C(t) denotes some continuous nonnegative function that 
does not depend on N. 

Next, we observe that, under Assumptions (Bl') - (B3'), the proof of Lemma 14,61 remains 
valid for time-dependent one-particle Hamiltonians. Thus, (|4.13p implies 

||M*) 1/2 9i(Wt)|| 2 < £(t) + (l + Mt)\\hnL~)(P(t) + 



Plugging this into (|4.33p yields 



Therefore, 



£{t) < 5(0) + jf ' ds C(s) (p(s) + £(s) + -IJj , (4.34) 

Next, we observe that, under Assumptions (Bl') - (B3'), the derivation of the estimate 
(I4,3ip in the proof of Theorem 14.11 remains valid for time-dependent one-particle Hamiltonians. 
Therefore, 

P(t) < 0(0) + jf ' ds C(s) (p( 8 ) + £(s) + . (4.35) 
Applying Gronwall's lemma to the sum of (|4.34|) and (|4.35p yields 

Plugging this back into (|4.35p yields 

which is the claim. □ 
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